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A FORMULA OF POLYNOMIAL INTERPOLATION.* 

By Webster G. Simon. 

Many proofs have been given of the following theorem of Weierstrass : f 
Let f{x) be a real function, defined and continuous for a ^ x ^b; and 

let e be an arbitrary positive constant. Then there exists a polynomial G(x) 

such that for a ^ x ^ b 

\fix) -Gix)\ < e. 

Of these proofs one of the simplest is due to Landau.J A proof will be 
proposed here, which, while closely related to Landau's, is perhaps more 
elementary, in that integrals are replaced by finite sums ; and at the same 
time it retains something of the elegance for which Landau's proof is 
conspicuous above others that are more elementary still. § The formula 
used here has some interest for its own sake, as a formula of interpolation, 
and as such may be compared, for example, with those of Borel,|| Faber,1T 
and Runge.** 

Without loss of generality, we may assume that < a < 6 < 1, 
since a linear transformation can be performed on x, if necessary, without 
affecting the essential conditions of the problem. Let the definition of 
f{x) be extended by setting f{x) = f(a) ior ^ x ^ a, f(x) = f{b) for 
6 ^ X ^ 1. Then f(x) is defined and continuous for ^ a; ^ 1. We 
shall now prove the 



* This paper, as originally written, contained in addition to the polynomial formula given 
here a trigonometric formula of approximation. While this paper was in the hands of the editors, 
Kryloff published in Bulletin des sciences mathematiques, second series, vol. 41 (1917), pp, 
309-320, a paper in which he establishes the validity of the same formula. So the reader is 
referred to Kryloff's article for the trigonometric formula. 

t Weierstrass, "tjber die analytische Darstellbarkeit sogennanter willkurlicher Functionen 
einer reellen Veranderlichen," Sitzungsberichte der Kgl. Preuss. Akademie der Wissenschaften, 
1885, pp. 633-639, 789-805. "Uber die analytische Darstellbarkeit sogennanter willkurlicher 
Functionen reeller Argumente," Werke, vol. 3 (1903). pp. 1-37. 

t Landau, "Uber die Approximation einer stetigen Funktion durch eine ganze rationale 
Funktion," Rendiconti del Circolo Matematico di Palermo, vol. 25 (1908), pp. 337-345. 

§ Cf., e. g., Lebesgue, "Sur ^approximation des fonetions," Bulletin des sciences mathe- 
matiques, series 2, vol. 22 (1898), pp. 278-287. 

II Borel, Legons sur les fonetions de variables reelles et les developpements en series de 
polynomes. Paris, 1905, pp. 79-82. 

It Faber, "Uber stets konvergente Interpolationsformeln," Jahresbericht der deutschen 
Mathematiker-Vereinigung, vol. 19 (1910), pp. 142-146. 

** Runge, Theorie und Praxis der Reihen, pp. 134-142. 
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Theorem. When n increases indefinitely, the polynomial in x {of 
degree 2n) 

t.m[i-{u-xYY 

(1) (Jn{x) = '— ;, , 

2 E (1 - tiY 

where ti = ijn, approaches the limit f{x) uniformly for a S x S b. 

First consider the case when f{x) is identically equal to unity, and 
write 

±[1- (i, - x)^]" 

(2) Tnix) = •^=5_,i , 

2 S (1 - t,T 

In this expression, the numerator is less than the denominator. For 
each of the two largest terms of the numerator is less tjjian or equal to 
unity, each of the two next in order of magnitude is less than or equal to 
(1 — ii^)", and so on, while in the denominator there will be a surplus of 
terms which will not be called into comparison in this way at all. There- 
fore 

(3) r„{x) < 1. 

On the other hand, each of the two largest terms of the numerator 
of (2) is at least equal to (1 — ii^)", each of the next two in order of magni- 
tude is at least equal to (1 — t^^)"-, and so on for at least k pairs of terms, 
if x is in the interval (a, b) and k, depending on n, is the greatest 
integer such that k/n is less than y, y being a positive quantity such 
that 0<a-y<b + y<l. Then 

2 i (1 - tiY 1 + i: (1 - t^y 



Tn{x) > 5 = 1 - 

2 E (1 - <i^)" Ed- tiY 

i=0 !=0 

For i> k,l - ti' < 1 - 7'. Therefore 

(4) Tn{x) > 1 S > 1 n • 

E (1 - U'Y Ed- ti^Y 

(=0 i=0 

Since evidently 0< 1 — t^< 1, n(l — r^)" approaches zero as a limit 
as n increases indefinitely. Now 

lim E (1 - m" = lim E (l - ^,)" = ^, 

as the following shows. 
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n 

Out of 2^ (1 — t^jn-y consider only the first [ V?i] terms.* Each of 
these terms will be decreased if n is put in place of i^, so that 

E(l - i'^ln'Y > [a/;^](1 - 1/n)" > K[^ln], 

1=0 

where K is independent of n, since (1 — 1/n)" approaches a limit, not 

n 

zero, as n increases indefinitely. Hence 2(1 — ii^)" increases indefinitely 

4=0 

with n. 

On combination of this result with (4) and (3), it is seen at once that 
to every e > there corresponds an N independent of x such that 

1 - e < Tn(x) < 1, n ^ N. 

Therefore when a ^ x ^ 6, t„(x) approaches unity uniformly as n in- 
creases indefinitely. 

Now let fix) be the arbitrary function for which the theorem is stated. 
Let 

E/WU - iu-xy]" 

(5) Vnix) = fix)Tnix) = — ^ . 

2 E (1 - tiY 

t=0 

It is an immediate consequence of the preceding that (pnix) converges 
uniformly to the limit f(x) . Thus, given a positive e, it is possible to find 
an Til, independent of x, such that 

(6) I (p„{x) — fix) I < e/3, n S tii. 

On the other hand, by subtracting (5) from (1) and taking absolute 
values, we see that 

1 1 fiti) - fix) I [1 - iU - x)^]" 

(7) I (r„(x) - cpnix) I ^ '-^ ^ . 

2 E (1 - UY 

<=0 

Since fix) is continuous on the closed interval S x ^ 1, we can find 
a 5 < 1 such that for any two values x', x", in the interval considered, 

\fix") - fix') I < «/3, when \ x" - x' \ < 8. 

It may be assumed further that 8 < y. The expression on the right of 
(7) can be broken up into two parts, as follows. Let the terms in the 
numerator for which — 8^ti — x^She added to form the numerator 

* [ Vrj] stands for the greatest integer in Vn. 
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of a new fraction S„' with same denominator, and let the remaining ones 
be similarly used to form a fraction Sn". It is assumed now that x is in 
the interval a S a; S 6. The sum S„' will be increased if e/3 is put in 
place of 1 /(if) — f{x) I . The resulting expression with the factor e/3 
removed contains only a part of the terms of t„(x), all of which are 
positive. Thus 

Sn' < e/3, 

in consequence of (3). In S„", [1 - {U - xYY < (1 - ^Y < 1, and 
replacing each difference | /(f<) - f{x) \ by 2M, where M is the maximum 
of the absolute value of /(a;), we obtain 

2M{n + 1)(1 - 62)n 



S„" < 



2 E (1 - k'Y 



The numerator approaches zero as n increases indefinitely, and the 
denominator becomes infinite. Then for a sufficiently large rh, inde- 
pendent of X, 

S„" < e/3, nmrh; S„' + S„" < 2«/3, n S n.., 
and 

I (Tn{x) — <Pn{x) \ < 2e/3, H g Ui. 

On combining this with (6), and taking Ho equal to the greater of ni 
and n2, we finally obtain 

1 <^n(a;) — f{x) \ < e, n S Ho. 

Therefore a-„(x) approaches the limit f(x) uniformly as n increases in- 
definitely. 

Chicago, Illinois. 



